Granular clustering self-consistent analysis for general coefficients of restitution 
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We study the equilibrium behavior of one-dimensional granular clusters and one-particle granular 
gases for a variety of velocity dependent coefficients of restitution r. We obtain equations describing 
of the long time behavior for the cluster's pressure, r.m.s. velocity and granular inter-spacing. We 
show that for extremely long times, clusters with velocity dependent coefficients of restitution are 
unstable and dissolve into homogeneous, quasi-elastic gases, but clusters with velocity independent 
r are permanent. This is in accordance with hydrodynamic studies pointing to the transient nature 
of density instabilities for granular gases with velocity dependent r. 
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I. INTRODUCTION 



Granular materials are present in many natural sys- 
tems and play an important role in our daily lives P, , 
and in the economy since an estimated 1.3% of the U.S. 
electric power consumption goes into grinding particles 
and ores Q . The interest in such systems ranges from the 
purely theoretical to dailypractical applications, such as 
in construction industry 0|. 

Typical granular systems (GS) are composed of large 
numbers of discrete macroscopic grains. Their shape is 
usually irregular but in the present work we will consider 
them as smooth regular spheres of diameter d in vac- 
uum, as a good approximation that still captures some 
of the essential physics of the problem. However, their 
dynamical and statistical properties may be affected by 
the presence of an interstitial medium, such as air or a 
liquid when their Bagnold number is small enough 0, • 
Granular materials behave in interesting ways exhibiting 
different features from ordinary solids, liquids and gases, 
such as arches redistributing loads to the sides of solid 
arrangements of grains and inelasticity-induced non triv- 
ial velocity distributions in rapid flowing granular gases. 
Even some basic laws of thermodynamics, such as the 
zeroth law, may fail when extended to GS pa, • 

A very important aspect of the behavior of GS is their 
inherent tendency to cluster, i.e., a compaction due to 
an enhancement of the rate of collisions inside the sys- 
tem, accompanied by granular cooling down (kinetic en- 
ergy reduction). Many authors have studied granular 
gas clustering, or compaction for denser systems, be- 
havior for GS which are initially homogeneous and start 
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with a given amount of kinetic energy 0, H, El El El i 

and subsequently left on their own to cool down. Their 
initial cooling behavior, the homogeneous cooling state 
(HCS), obeys Haff's Law [H IT^ for the granular tem- 
perature (the typical internal average kinetic energy for 
the grains). 

For longer times and rarefied granular systems, the 
inter-granular collisions tend to correlate the motion (ve- 
locities and positions) of the grains and techniques from 
Kinetic Theory will not be reliable anymore, at least in 
its simplest form 0, 0, El . Theoretical models have 
developed upon this notion and obtained scaling forms 
for the transport coefficients 0,0] when the velocities 
of masses of grains become correlated inside a region of 
a certain characteristic length. 

For clusters coalesced from smooth granular gases (no 
tangential restitution) there are no mechanisms for the 
exchange of angular momentum and rolling does not oc- 
cur (a totally irrelevant sliding motion between grains 
may still be present, but the rotational motion of the 
grains is not coupled to the translational motion if the 
grains are smooth). It is legitimate, from a theoretical 
point of view, to ask whether such structures are really 
permanently stable or some other mechanism could lead 
the system again to gaseous homogeneity. 

It has been known, for quite some time, that the hydro- 
dynamic approximation predicts density instabilities for 
inelastic, smooth, hard-sphere granular systems at zero 
gravity [18j (velocity independent coefficients of restitu- 
tion) while for equivalent systems with velocity depen- 
dent coefficient of restitution (r) , such as the viscoelastic 
model ,19]] the instabilities are only transient j2jj. Thus, 
the dependence of the coefficient of restitution on the ve- 
locity might be the cause of a possible cluster break-up. 
That is due to the fact that r tends to 1 as the impact rel- 
ative velocity tends to zero. However, to simulate a clus- 
ter break-up can become computationally very costly, if 
no approximations are used. In most event driven Molec- 
ular Dynamics simulations the coefficient of restitution 
has to be set to 1 as the relative velocity becomes smaller 
than an elastic threshold, so that collisions at relative ve- 
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locities lower than the threshold will be elastic [Til |20| . 
In order to be more realistic, one has to be able to reach 
extremely long times and small velocities. This is the 
goal of the present work. We will study a simple qual- 
itative model (a 1-dimcnsional cluster) and obtain the 
asymptotic solutions for its behavior at extremely long 
times, not yet accessible to computer simulations. 

We study systems with coefficients of restitution given 
by the general form [2l| , at small relative velocity g 



A | — 

9o 



(1) 



where < A < 1, and the constant go sets the dissipation 
scale, and m > 0. The cases m = (velocity independent 
r) and m = 1/5 (viscoelastic) are particular cases of 
Eq.H 

We assume that the system starts as a granular gas 
inside a bounded 1-dimensional container by walls at the 
extremities: one elastic and one inelastic. A cluster forms 
at the inelastic end and a gas still exists at the elastic 
side after some initial large time interval 0, ■ The 
granular gas portion will consist ultimately of a single 
grain. We will use this assumption in what follows but 
this is not really essential, and our results could easily be 
extended to multi-grain gases, provided the gas kinetic 
energy goes asymptotically to zero. 

We will obtain equations of motion for the dynamic 
and quasi-thermodynamic granular quantities: gas ve- 
locity, the cluster granular temperature (in fact its r.m.s. 
velocity) and internal pressure. We will show that the 
cluster is unstable at extremely long times, except for 
the case of velocity- independent case (m = 0). Purely 
dynamic effects, due to energy dissipation, can be respon- 
sible for cluster formation but can not hold it together in- 
definitely, except for the unrealistic velocity-independent 
coefficient of restitution case. 

This paper is organized as follow. In Section II, we 
write the basic equations of collision dynamics. In Sec- 
tion III, we present our model and its corresponding 
equations. In Section IV, we study the stability condi- 
tions for the cluster. In Section V, the long time behavior 
is extensively analyzed. In Sections VI and VII, we dis- 
cuss some of the consequences of the model and make our 
concluding remarks. The more technical aspects of this 
work are to be found in Appendices IVIIII and IIXI 



II. BASIC ASSUMPTIONS 

Our approach is qualitative in nature and we do not 
expect it to reproduce the detailed behavior of a true 
granular cluster but only a few of its crudest properties. 

Our model consists of an unforced one dimensional sys- 
tem with N + 1 identical smooth grains, with unit mass 
M = 1, which can be modeled by inelastic hard-spheres 
with velocity dependent coefficient of restitution r, given 
byEq.n 



The grains are confined in a region of length L + (N + 
l)d (d is the granular diameter) by an elastic wall on the 
right and an inelastic one on the left, as shown in fig- 
ure n Relative to each other, the elastic wall presents 
some similarities to a "hot" wall while the inelastic one 
would be a "cold" one, as in reference |23, since the in- 
elastic one is where the energy is taken out of the system. 
However, we are not injecting any amount of energy into 
the system and a steady-state does not develop. The in- 
elastic collisions at the inelastic wall are governed by the 
same Eq. ^ the wall being an infinitely heavy grain. The 
velocities of two colliding grains will be (V{, V^') after the 
collision, and (Vi, V2) before it. They are related by: 



v; = 



V' 



l-r 



Vl 



1 + r 



Vn 



(2) 



(3) 



The system starts with a given initial amount of kinetic 
energy that will be dissipated due to the internal colli- 
sions. A partial clustering of grains will initially occur 
at the inelastic wall due to the pressure of the remaining 
gas [22^. The cluster phase is formed when the relative 
velocities are large compared with g and the coefficient 
of restitution differs appreciably from 1. The gas will 
lose particles to the cluster until only one gas particle 
is left 15]. The difference between clustering and col- 
lapse is important, since in the later an infinite number 
of collisions occur among the particles in a finite amount 
of time. However, the cluster is only a very dense con- 
centration of grains with small relative speeds. In one 
dimension clustering can precede collapse (for r close to 
1) 0. 

After some initial interval of time, the velocities of all 
grains, cluster's one and the gas one, will be much smaller 
than the inelastic velocity scale go in Eq. ^ 

We label the gas particle as the th particle while the 
cluster ones are labeled from 1 to N. We assume go ^> 
vo 3> Ui=i ... } n- Although the velocity of the gas particle 
is much larger than the velocities of the particles forming 
the cluster, the scaling factor g will also be a lot larger 
than the velocity of the gas particle for long times. The 
following hierarchy for the small expansion parameters 
holds (for m > 0): 



» 



(4) 



The logic of Eq. 01 is that at sufficient long times 
— — > and (since the gas will keep pumping energy 
into the cluster) the ratio Vi /vo will not vanish (that will 
be checked a posteriori, with respect to the long time be- 
havior of the cluster). We can assume (for m > 0) that 
K/t>o| ^ |^o/go| m - The second inequality comes from a 
choice that the time we choose to start our calculations 
is large but not too large, since for even larger times, 
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certainly \vi/vo\ 3> |t>o/go| m - The last inequality is a 
direct consequence of this choice of initial time. Higher 
order terms shall be discarded. 

We assume that our model has initial conditions satis- 
fying Eq.H 



III. GAS-CLUSTER EQUILIBRIUM 

A. Gas pressure 

In order to calculate the pressure exerted by the gas on 
the cluster we need to take into account the momentum 
exchanged between the gas particle and the cluster af- 
ter each collision between gas and cluster. A gas-cluster 
collision is completed when the gas particle leaves the 
cluster with an outgoing velocity (which is close in ab- 
solute value to its velocity before collision) much larger 
than the typical cluster's particle's velocity. This is il- 
lustrated in Fig. H We notice that as the gas particle 
collides with the cluster, it is as if the "fast particle" 
(gas) would pass through the cluster's ones, collide with 
the inelastic wall and go all the way back in the direction 
of the elastic wall. The process can be repeated until the 
"fast particle" crosses the whole cluster, back and forth. 
Thus, we obtain, to leading order, the momentum ex- 
changed between the gas and cluster. This is calculated 
in Appendix IVIIII Eq. 0H1 (notice that V is the speed of 
the incoming gas particle, V > always): 



V = -(2iV + l)A 



V 2 , 



(5) 



where the nipping of the "fast particle's" velocity at the 
inelastic wall is correctly taken into account in the final 
result above. 

From Eq. ^2(m > 0), we obtain 



AR 



Pgas 



cluster 



At 



(l + m)V . 
2 V £ ' 



(6) 



where p gas is the pressure exerted by the gas on the clus- 
ter, and AP c i ust er is the inelastic change of cluster mo- 
mentum due to the collision with the gas particle. 
For m = 0, Eg. EB1 gives us 



Pga 



l-(l-A) 

1 + (1- A) N 



V. 



(7) 



We will assume in the following that the product NA 
is small. Our results, being valid for small NA in the 
to = case, will show that collapse happens in a finite 
amount of time. This will be certainly valid for the case 
of large NA. 



B. Cluster's variables 

Other important variables exist that we need to take 
into account. The variance of the cluster's particle's ve- 



locities a 2 is one of them. It is defined as 

.. -.2 



N 



(8) 



where the cluster's center of mass velocity is given by 



N 



(9) 



Another such variable is the mean granular spacing, e. 
It is defined by (with xn+i = 0) 



2^=1 \ x i ~ x i-i - d\ 



N 



(10) 



The cluster's full size is thus N(d+e). A somewhat crude, 
but reasonable, approximation for the cluster's center of 
mass velocity can be obtained by assuming that the clus- 
ter expands or contracts rather uniformly (at least for 
small e), and thus 



JV 



Nii = 



N(N +1) 

27V ' 



(11) 



In the limit of large N the expression above reduces to 

N 



The center of mass' acceleration then reads 
N(N + 1) 



2N 



-e. 



(12) 



C. Wall pressure and mean spacing 

The cluster feels two external sources of pressure: the 
pressure from the gas and the pressure due to interactions 
with the inelastic wall. That pressure can be calculated 
using a crude approximation by assuming that the mo- 
menta exchanged between the particle labeled N at every 
collision is of the order of 2a and the rate of collision is 
a/e [25| . The mean-field wall pressure is then 



2a 2 

Pwall = > 0. 



(13) 



The equation of time evolution for the mean spacing e 
is obtained by writing Newton's law for the cluster (p wa u 
is a positive force while p gas is a negative one): 

^cm — Pwall ~~\~Pgas- 

Using Eas. Iffl Upland 1131 we obtain the mean- field equa- 
tion for to > 0: 



N 2a 2 (l + m)V 

—e = h '-— e. 

2 e 2 V 



(14) 



The equivalent equation for the case to — is obtained 
from Pg as given by Eq. 1461 

2 £ - e + 1 + (1-A)» V - (15) 
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D. Energy dissipation inside the cluster 



IV. DIMENSIONLESS ANALYSIS 



A calculation similar to that for the gas velocity re- 
duction is done in Appendix IIXI for the energy of the 
cluster, due to the effect of gas-cluster collision. At the 
order of approximation we have set, |V/<7o|"\ the correc- 
tions will be of the order & 2 \V/go\ 2m (much smaller than 
cr 2 /e). The comparison of the cluster's kinetic energy be- 
fore and after a collision with the gas particle then reads 



JV 



/V 



£K"') 2 = £K) 2 



(16) 



The cluster's kinetic energy is not affected by the gas- 
clustcr collision, at the order of approximation used. This 
tells us that only internal collisions will be important for 
cooling down the cluster. 

In a mean field approximation, the energy dissipated 
per particle corresponds to the product of the energy lost 
in each collision and the rate of collision per particle. 
For a A^-particle one-dimensional gas confined in a free 
volume 1\d, with the typical velocity variance a (v is 
the typical velocity, of the order of the square root of 
the variance a), the energy loss per collision corresponds 
roughly, in dimensionless terms, to 



Av oc 



il+m 



The rate of collision is proportional to v/1%d |25|- So, 
in order to calculate the rate of dissipation of energy in- 
side the cluster, with typical inter-spacing becoming very 
small as t — > oo, but not yet zero, we need to estimate the 
internal collision rate and multiply it by the loss of en- 
ergy for each collision. The rate of collisions for a cluster 
only differs from that of a gas because its inter-spacing e 
is very small compared to L. Taking it into account, the 
cluster's collision rate is now obtained as 



q = N 



2s 



The variation of the cluster's typical internal velocity is 
due to the collision between grains. We assume that col- 
liding grains will have relative velocities of the order a 
and the quasi-elastic collisions (in the limit g — > 0) will 
switch those velocities with a small loss 



Act = -A — 
\9q 

The change in a per unit time is the product qAa. Thus 
we obtain the heuristic equation: 



NAcr 2 
2s~ 



(17) 



The equivalent form for m — is derived from Eq. 1531 
from appendix B: 



We want to study, qualitatively, the conditions for the 
cluster to be stable and shall use the approximate equa- 
tions of motion obtained above. However, we will not be 
interested in the fine details of the equations themselves, 
only in their asymptotic behavior in time. Then, we will 
rewrite Eas. 151 PHI 1151 and 1171 in a completely dimension- 
less form as (notice that V,a, e > for all t). First, the 
equations for V and a have the same form for both m = 
and m > 0. They read 



V = -V 2+m , 

a 2+m 
a = . 



The equation for e for m = reads 

a 2 ■ 
i = — +V, 

e 



and that for m > reads 



V 



(19) 
(20) 

(21) 

(22) 



In order to recover the dimensional units, remember 
that a and V are given in terms of go, and e is measured 
in terms of L/N. A few dimensional constants have to 
be used in Eas. 11912*2*1 in order to make both sides dimen- 
sionally coherent. 

Eq. [""J"] can be exactly solved, obtaining 



V 



V 



(l + (l + m)V a 1+m t) 



(23) 



This is the extension of Haff 's law 0, |2|| to the cases 
described by Eq. [2 We observe that as 



V 



For the case of velocity independent coefficient of restitu- 



tion, m = 0, T q 



f- 



For the viscoelastic coefficient of 



restitution case pjj, m — 1/5, T g ~ t~3 as expected [26j . 
Eqs."3Ul l22land"2*Tlwill constitute the system to be solved 
in the following. 



V. LONG-TIME BEHAVIOR 

We must keep in mind that there is an implicit velocity 
scale go that divides the velocities variables whenever a 
power of m comes into play (a consequence of the form 
of the coefficient of restitution) . For the constant coeffi- 
cient of restitution case, we obtain a useful equation by 
combining Eq.E3]with Ea. [2*1*1 



(18) 



V 



c - 



(24) 
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where the constant cq is related to the initial conditions 
for V, a and e by: 



c = £o - Vq + 



1 



(25) 



The equations describing the granular cluster are valid 
in the limit when V 3> cr,i. In the following we analyze 
the cluster behavior for the whole range of values of m 
based on Eas. 11912*41 The figures are obtained by solving 
numerically the Eqs. 1191241 



A. m = 

In this case, the coefficient of restitution does not de- 
pend on the impact relative velocity. We have: 



co = £q — Vb + o" 
i = V-a-V . 



-V 



(26) 
(27) 



As t — > oo, we observe that V and a tend to zero |27| 
and thus e f» — Vo < 0. 

In practice, it is impossible to observe the asymptotic 
limit above since the collapse happens on a finite amount 
of time (see figure [2J|. The cluster is thus stable and will 
not dissolve itself. 



B. m > 

When the coefficient of restitution depends on the ini- 
tial relative velocity, i.e. m > 0, we can see that the 
physical behavior of the system changes qualitatively, as 
shown in Appendix IVIIII 

A scaling argument can be used and compared with the 
result of simulations in order to obtain the very long time 
behavior of the variables V, a and e. The asymptotic 
solutions for s and a can be written as powers of time 
and log-time: 



V-t? 1 , cr ~ t 132 (In t) a \ 

where we have already determined (3\: 

1 



^(lnt)" 3 



(28) 



01 = -. 



1 + m 

The solutions for 02, &3 and 03 are 

1 

«2 = , 



(29) 

til 

02 = 0, (30) 

a 3 = - — , (31) 
m 

03 = 1. (32) 
Thus, the long time behavior of a and s is given by 

cr~(hii)~™, and e ~ i(lni) - ™. (33) 



These are self-consistent, logarithmically corrected so- 
lutions for a and e at long times. However, we need to 
look into the long-times behavior of a with more detail. 

For m > there is no granular collapse. After some 
transient time, the cluster will grow almost linearly (as 
can be seen in figure"""] consistent with the main behavior 
of e ~ t) and will eventually occupy the whole container, 
in fact becoming once again a granular gas with inter- 
spacing e ~ L/N. It takes an enormous amount of time 
for this to happen. This is illustrated in figure |3| where 
we compare the cases m = 0.2, m = 1 and m = 2. 

At long-times, our model becomes quasi-elastic, to a 
very good approximation, when m > 0. For the velocity 
dependent case, the internal dissipation for the cluster 
becomes negligible (02 = 0) but its internal energy is not 
conserved. The apparent contradiction between 02 = 
and (7 — > does not hold since a decays as a power of 
lni. Even more significantly, in our model e < L/N and 
the growth of e has to be cut-off correspondingly. Since 
our mean-field equations do not impose a boundary to 
e, then Eq. 1201 will give us a — > as t — » 00, consistent 
with /?2 = 0. In reality, after reaching the cut-off size, 
normal gas dissipation takes over and the former cluster 
will follow Haff 's law for energy dissipation again. 

Another important consistency argument can be ex- 
tracted from Eq. If we take the limit m — > before 
the limit t — > 00 is taken, we observe that e = a = re- 
sults. This is in complete agreement (remember that the 
initial time is taken to be long for the m > case) with 
our result that a granular collapse happens in a finite 
interval of time when m = 0. 

It is interesting to notice at this point that in refer- 
ence [2(| the authors simulate a velocity dependent gran- 
ular system with an elastic threshold (r = 1 below a 
certain threshold relative velocity) and supposed their 
results to be extensible to the viscoelastic regime. This 
is in accordance with our results. However . the form of 
the coefficient of restitution in reference [2(| much more 
closely mimics the m > 1 than the case m — 0.2. As 
we observe in our calculations, e will neither tend to zero 
nor remain stable in both cases, which is consistent with 
the results in reference 201. 



C. Oscillations 

An interesting feature we have observed are very low 
frequency size-oscillations, at very long times. Eq. QUI 
predicts oscillations with decreasing frequency. We can 
study the case of a small perturbation on e such as: 

e = t(]nty™ [1 + 0] . 

We obtain, after some straightforward algebra, the 
asymptotic equation for the relative perturbation <j>: 



hxt 



2 hit 



t 2 



-cf> = 0. 



(34) 
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It is similar to a low frequency damped harmonic oscil- 
lator, with a frequency that goes to zero as i _1 (ln£)2 . 

The effect of small perturbations in the asymptotic 
value of e is rather hard to observe directly. However, 
we observed it by initially running our simulations in or- 
der to obtain aged values of e, a and its derivatives. We 
then perturb e as (1 + A)e, with A = 1.0 x 1(T 4 . 

We run two subsequent calculations, with an aged and 
unperturbed solution as the initial condition, and an- 
other for the perturbed one. Their difference should also 
obey Eq. |^ The result is plotted in figure U in a loga- 
rithmic scale (we plot the absolute value of <f>; the signs 
correspond to whether <j> is positive or negative). 

It can be seen that the period is indeed increasing (it 
is of the order of the total time, consistent with a "fre- 
quency" of order i _1 (lni)2. 



VI. CONSEQUENCES OF THE MODEL 

The most immediate consequence of the present model 
(for to 7^ 0) is the evidence it provides of the transient 
nature for some of the granular singularities in a freely 
cooling granular gas with velocity-dependent coefficient 
of restitution. This indicates that a hydrodynamic treat- 
ment might be adequate for such systems, at least after 
a transient time. Also, we deduce from our results that 
purely dynamical effects cannot give rise to permanent 
clusters if m ^ at the zero-energy feeding regime. 

Another consequence is the eventual evaporation of 
clusters for smooth granular systems. The inviscid Burg- 
ers' equation has been proposed as a mechanism of for- 
mation for a granular cluster with velocity independent 
coefficient of restitution [Tlll2S| . For systems with m/0 
one may ask whether that equation is still adequate, and 
what kind of regime might replace it, in the evaporative 
period (at extremely long times). Work is currently un- 
der way along this direction. 

The non-collapse when m ^ gives us hope that 
it might be possible to treat two- or three-dimensional 
clusters as a very dense, but non-singular, granular 
phases (for smooth systems) describable by internal, non- 
diverging, variables (maybe even similar ones to the a 
and e used in this manuscript). That could make it easier 
to incorporate the treatment of clusters into the hydro- 
dynamic methods available today. 



VII. CONCLUSIONS 

We study the long term stability of unforced granu- 
lar systems, in which clusters form, with the help of a 
qualitative, microscopic model that makes it possible to 
look at clusters at extremely long times, not available to 
computer simulations. 

We assume a general form for the coefficient of resti- 
tution that includes the well known velocity-independent 
and viscoelastic models as special cases. 



We are interested in this problem for two main rea- 
sons. Firstly, despite its apparent simplicity, a granu- 
lar cluster's behavior, at extremely long times, depends 
on the amount of inelasticity (which can be defined as 
q = ^-)- According to our model, if the coefficient of 
restitution becomes 1 as the relative velocity of impact 
tends to zero, as with most realistic systems, then clus- 
ters of rigid, smooth spheres will be unstable (at least at 
zero gravity) . This suggests a rich dynamical behavior for 
our granular gas that comprises an initial homogeneous 
phase in which Haff's law predicts the evolution of 
the average granular temperature. The system goes into 
phase separation after a transient time and the global 
kinetic energy varies with a different power of time • 
After a very long waiting time, the external granular gas 
pressure no longer keeps the cluster particles together 
and the cluster finally dissolves into an extremely slow 
moving homogeneous granular gas. Haff's law will once 
again apply to this gas (since to 0). This is not in 
contradiction with the results in reference [2g since the 
results therein apply to systems with velocity indepen- 
dent coefficients of restitution (to = 0) . 

Secondly, for velocity dependent coefficients of resti- 
tution, the clusters are not truly collapsed, but behave 
instead as very dense, fluid phases (for zero surface fric- 
tion and zero gravity). In fact, we could think of the 
gas-cluster phase coexistence boundary as a smooth sepa- 
ration between the granular gas and cluster phases, with- 
out a singular boundary, except for the case of constant 
coefficients of restitution (to = 0). An appropriate con- 
tinuous hydrodynamic treatment for it might be possible. 

Questions arise concerning the long times dissolution 
of granular clusters: will they obey the same equations 
as the ones that are found to apply for the collapsing 
phase? Since the irreversibility of the "microscopic" , e.g. 
granular, dynamics prevents time-reversal to apply, the 
dissolution equations might be quite different from the 
collapse ones. This is yet to be understood. 



VIII. GAS-CLUSTER MOMENTUM 
EXCHANGE 

A. m > 

From Eqs. [21 and |3| we obtain the post-collision veloci- 
ties for two particles of the same mass (the collision time 
is taken to be zero): 



A 
2" 



vi - v 



9o 



A ui — V[) 

V + — ; Vi, (35) 



Vi - Vq 



9o 



A 


V% - V 


( A 
vo+[l- ¥ 


vi - v 


rn \ 








J Vi 


~2 


90 




go 





vi. (36) 



The final velocity of the "fast particle" after crossing 
the cluster can be calculated by the equations of basic 
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collision dynamics. Assuming go 3> \vq\ =F» |ui|, we 
can expand the last term of Eq. Q]as follows: 



9o 



vo 


m 


i - 




m 


9o 




uo 




V 


m 










- m 




9o 











(37) 



The velocities of particle 1 and the gas particle, after 
the first collision, can be rewritten with the help of Eg. 1551 
and Eg. 131)1 as: 



where p" cl = J2Zi v i- 

The momentum received by the cluster due to these 
collisions is then 



^■Pcluster2 



= N- 



N- 



4 



V' - (1 + m) 



A 



Pel, 



V - (1 + m)- 



.4 



P% 



Before the collision with the elastic wall, the fast par- 
ticle velocity is given by: 



-V 



A 
~2 



V+{l + m) 



A 



vi, 



now being the fast particle, and the gas one is now slow 
(primes stand for fast-slow collisions): 



A 
~2 



V - (1 + m) 



.4 



Vi. 



After i collisions, the fast particle velocity will be the 
l-Xh one 



-V + l 



A 



V + (1 + m) 



A 



i.=i 



and the [I — l)-th particle (which suffered two collisions) 
has the velocity: 



A 
~2 



V - (1 + m)- 



.4 



After colliding N times, the fast particle will reach the 
inelastic wall. Its velocity, prior to that collision, will 
read then: 



v' N = -V + N 



A 



V + (1 + m) : 



.4 



Pel, 



A 

v = V -N- 



V" + (l + m)^ 



Pch 



V-(2N+1)~ 



+(l + m) 



.4 



V + 

(Pel - Pel), 



After that last collision, the gas particle has the veloc- 
ity: 



= -y + (2iv + i)4 



-(1 + m)- 



V - 

(Pel-Pel), 



Hence, the absolute value of the gas particle velocity 
varies (for a single gas-cluster collision cycle) as: 
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-(2JV + 1 



.4 



+(l + m)- 
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(Pel - Pel), 

71 

V, 



(38) 



where p c i = 53fcLi v i, the cluster's total momentum be- 
fore colliding with the gas. 

Hence the total momentum given by the gas to the 
cluster (first part) 



Ap, 



lu sic 



-N- 



V - {1+m 



A 



Pel- 



After the collision with the inelastic wall, it reads: 



V = v'j, 



V-(N + 1)^ 



A 

V-(l + m)- 



Peh 



The procedure for the calculations of how the fast particle 
traverses the cluster is similar to the one above and the 
final result is (up to the same approximation order) 



= V'-N- 



V' + (l + m)- 



P"el, 



where we discarded terms of the order C(| V/go| "*)■ 

Thus, the total momentum absorbed by the cluster 
from the gas-cluster collision is given by 



Ap, 



cl tot 



^Pclusteri ^Pcluster2 

v 



-(1+m)- 



= -(l + m)A 



90 
V_ 
9o 



(Pel + Pel) 



Pel- 



(39) 



The reader should notice that we assume the initial time 
to be large enough so that quantities such as NA\V/go\ m 
are small and the total dissipation per gas-cluster colli- 
sion can be a small fraction of the gas kinetic energy. 

Equations 1551 and |5J5]are the fundamental result of this 
appendix. We can transform them into rate equations by 
determining the rate of gas-cluster collisions. The time 
interval between successive collisions is given by At = 
2L/V. 
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We obtain the equation governing the behavior of the 
absolute value of the gas velocity: 



V 



(40) 



The equation for the gas pressure, the rate of transfer 
of momentum is also obtained 



Ap 



'cl tot 



At 



- (1 + m) 2L 



V Pc 



Notice that for a large cluster, N 3> 1, we can write 

Q- + m)V • IMS 
Pgas = „ ' ~ e. (41) 



V 



The gas pressure in this case will be given by (similarly 
to case for m > 0) 



Pgas 



a Pc _ ((i-A) N -iy 

At 2L 



(45) 



We can see that the gas pressure is related to the 
change in gas velocity through 



_ l-(l-A) 
Pgas - ! + (1 _ A ) 



A 



JV 



V. 



(46) 



There is a clear change in the gas pressure regime for 
m > compared with the more commonly used case 
of m = 0. This makes the pressure applied by the gas 
weaker since, for m > 0, the factor V (see Eq. 146(1 is 
multiplied by a factor i/V (see Ea. l4T|) . 



B. m = 

The case of a constant coefficient of restitution deserves 
a separate treatment. In this case we assume that A <C 1 
and r = l- A»l. 

After a collision with a slow particle, the fast particle 
acquires a velocity v' = (1— A)v. At the end of a sequence 
of N such collisions, the velocity of the fast particle (be- 
fore colliding with the inelastic wall) will be 

v N = (l-A) N v . 

The momentum exchanged with the cluster is then 

A Pcl = ((l-A) N -l)v . 

After colliding with the inelastic wall, the gas particle 
has a velocity vn = —(l — A) vo- After colliding another 
N times with cluster particles, the gas particle velocity 
will be 

v'^ = -(l-A) 2N v a . 

The momentum exchanged with the cluster this time is 
then 

Ap c2 = -(1-A) N ((l-A) N -l) v . 

The total change in velocity for the gas particle, after 
collision with the elastic wall, is given by 



AV 



(l-(l-Af N )V, 



(42) 



The rate of change of V is given (see the calculation for 
m > above) 



V 



1-0_A) 
2L 



2 A 



V 



(43) 



The total momentum gained by the cluster after the 
collision is then: 



Ap c = Ap cl + Ap c 



((1 - A) N — l) V. (44) 



IX. CLUSTER ENERGY DISSIPATION IN THE 
GAS-CLUSTER COLLISION 

A. m > 

In order to show that the cluster's kinetic energy is 
not affected by the gas-cluster collision on the order of 
approximation we have chosen, let's consider the sum of 
the velocities after the first passage of the gas particle all 
the way to the inelastic wall: 



N-l 
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m- 
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90 
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NV. 



(47) 

Our calculations are carried out to order N\V/go\ m in 
the development of the coefficient of restitution, in the 
spirit of Eq. We suppressed terms coming from orders 
smaller than N\V/g \ m - 

Let's also consider the sum of the velocities after the 
passage back of the gas particle 
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£' 

i = l 



1 - (1 + m) 
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JV— 1 



A 



i=0 



These equations can be added up giving 



NV. 

(48) 



A 



5>r-5> = -(i+mM - ( 49 ) 



JV 



yielding the pressure exerted on the cluster by the gas. 

We can proceed along similar lines for the sum of 
square velocities and obtain 



JV-l 



JV 



E«) 2 = E^ 2 -^- E«" 



i=0 



and 



A 



JV-l 



V 
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E«") 2 = E«) 2 + AV t E 



i=0 



i=l 



m JV-l 



(50) 



(51) 



i=0 
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Equations El HOI an d El show that the kinetic energy 
of the cluster remains the same: 

N N 

£(O 2 =£^ + 0(" 2 Wffo| 2m ). (52) 

i=l i=l 

B. m = 

As shown in appcndix lVIII Bl the gas grain pumps mo- 
mentum into the cluster. We will assume that NA <C 1. 
This is not too restrictive to our argument since we will 
show that a long-time granular collapse happens for the 
quasi-elastic velocity-independent coefficient of restitu- 
tion case. Thus, it will certainly happen for the case 
when NA is large too. 

We noticed that as the fast grain collides with the clus- 
ter it gives energy to it by changing the particles' veloci- 
ties by an amount proportional to NAV. After squaring 
all cluster particles' velocities (relative to the center of 
mass of the cluster) , adding them all up and subtracting 
the initial value of it, we obtain a rate of energy, pumped 



into the cluster, proportional to the product of V and i. 

The rate of change of a has two main contributions: a 
negative one from internal collisions; and a positive one 
from gas-cluster collision. The second one is negligible 
and we do not take it into further account. The reason 
for it goes as follows. Since \V\ \iV\, if we assume 
\iV\ > a 1 je then a will decay more slowly when the 
energy pumping term is present but the wall pressure 
term in Eq. 1211 will still be much smaller than the gas 
one. However, we can check a posteriori that even when 
the energy pumping term is not present, the mean inter- 
spacing falls at a linear rate in a finite time (collapse, see 
figure [21 • It yields 

a 2 /e~ \V\ > \Ve\. 

Thus, we only keep the internal collisions dissipation 
term in the equation for &: 

2 

<rw- — . (53) 
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FIG. 1: The iV particles of the cluster and the gas particle. 
The average distance between two consecutive particles of the 
cluster is e (not shown), the diameter of each particle is d, and 
the total length is L + (N + l)d. 
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FIG. 3: Time evolution for a few examples of velocity depen- 
dent coefficients of restitution: full line m = 0.2 (viscoelastic 
model); dashed line m = 1; dotted line m = 2. 



FIG. 2: Coefficient of restitution independent of the velocities. 
The only stable case at all times. The mean-spacing e is 
measured in units of L/N. 
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FIG. 4: Plot of the unperturbed value of e minus the per- 
turbed one, as a function of time, normalized by e itself. It 
is consistent with Eq. I34I a damped-harmonic-oscillator-like 
equation. The sign corresponds to whether the oscillation has 
a positive or a negative value, since the plot corresponds to 
the logarithm of the absolute value of the difference between 
the two freely evolving solutions. 



